Enhanced Thermoelectric Figure of Merit in 
Edge Disordered Zigzag Graphene Nanoribbons 



H. Sevingli^ and G. Cuniberti^ 

^Institute for Materials Science and Max Bergmann Center of Biomaterials, 
Dresden University of Technology, D-01062 Dresden, Germany. 

We investigate electron and phonon transport through edge disordered zigzag graphene nanorib- 
bons based on the same methodological tool of nonequilibrium Green functions. We show that edge 
disorder dramatically reduces phonon thermal transport while being only weakly detrimental to 
electronic conduction. The behavior of the electronic and phononic elastic mean free paths points 
to the possibility of realizing an electron-crystal coexisting with a phonon-glass. The calculated 
thermoelectric figure of merit {ZT) values qualify zigzag graphene nanoribbons as a very promising 
material for thermoelectric applications. 
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The isolation of graphene [1], a one-atom thick sp^- 
bonded planar carbon sheet, and quasi one-dimensional 
graphene nanoribbons (GNRs) have inspired research in 
many directions [2] . GNRs are semiconductors with band 
gaps depending on their width [3, 4]. The magnetic as 
well as electronic properties of GNRs are strongly depen- 
dent on edge shapes and very sensitive to width variations 
throughout the ribbon [5]. GNRs with average widths 
less than 10 nm are fabricated [6], and control of GNR 
edges with an atomic precision has not been achieved yet, 
so edge disorder is an intrinsic property of GNRs. 
Charge transport through edge disordered GNRs has 
been studied in detail in the last few years [7, 8] . In par- 
ticular, Areshkin et al. showed that GNRs with armchair 
edges are extremely sensitive to edge disorder whereas 
for those with zigzag edges (ZGNRs), the outstanding 
transport properties are weakly affected within the first 
conductance plateau (FCP), i.e. the energy spectrum be- 
tween the second conduction and valence band edges [8] . 
For thermoelectric energy conversion, on the other hand, 
a low thermal conductance is required and indeed high 
thermopower values up to 100 /iV/K is reported for 
graphene [9] . However two-dimensional graphene has ex- 
tremely high thermal conductivity which is dominated 
by phonons [10], which nevertheless is strongly affected 
by graphene edges [11]. The question is whether it is 
possible to overcome the high thermal conductivity so 
that GNRs might be good candidates for thermoelectric- 
ity. Recently, significant reductions of phonon transport 
through nanotubes due to isotopic or Anderson-type dis- 
order were reported [12, 13, 14, 1.5]. It is also shown that 
phonon thermal conductance through Si-nanowires can 
be reduced by up to two orders of magnitude due to sur- 
face roughness or surface decoration which give rise to 
high thermoelectric coefficients [IG], as it was predicted 
theoretically that the thermoelectric figure of merit can 
be enhanced in low dimensional systems [17]. 
In this Letter, we investigate electronic and phononic 
transport properties of edge disordered ZGNRs on an 
equal footing. We show that Kph can be reduced dra- 



matically while the electronic transport can stay almost 
intact at the FCP, making edge disordered ZGNR an 
electron-crystal at these energies whereas a glass transi- 
tion for phonons is achievable. 

Methods— We identify ZGNRs with the number of 
zigzag chains, Nz, which also determines their widths. 



(a) 



h L- 




H 




^^^^ 


^^^^^^^^^^^^^^^ 



LEFT 



CENTER 



RIGHT 



(b) 



(c) 




1600 



400 800 1200 

Energy (cm ') 



1600 



FIG. 1: (Color online) Schematics of the system (a). The 
central region with edge disorder has length L, and is con- 
nected to semi-infinite ZGNRs free of disorder. Transmission 
spectra of phonons, Tph, for pristine and disordered ZGNR 
are plotted for Nz —10 and 20 in (b) and (c), respectively. 
The dark (red) regions represent the difference in transmis- 
sion between pristine ZGNR(IO) and ZGNR(8) (b), and the 
same for ZGNR(20) and ZGNR(18) (c). Tph is suppressed 
with increasing L =15.7, 63, 252 and 504 nm. Elastic mean 
free paths are given in the insets. 
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We use the common partitioning scheme for both elec- 
trons and phonons by dividing the system into three re- 
gions, namely left, right and central regions [18, 19]. The 
central region includes the disordered part, whereas the 
left and right regions are taken as semi-infinite perfect 
GNRs of given width (Fig. 1). 

For constructing the dynamical matrix, we use the fourth 
nearest neighbor force constant approximation (4NNFC), 
which yields phonon dispersions in agreement with den- 
sity functional theory (DFT) calculations for graphcne 
and carbon nanotubes [15, 20, 21]. For the case of GNRs, 
modification of force constants for the edge carbon atoms 
will improve the results [22, 2.'!], in the sense that it will 
result in a blue shift in the phonon density of states. 
Therefore a subsequent reduction of lattice thermal con- 
ductance and a further but minor improvement to our 
results in reducing lattice thermal conductivity can be 
expected, but it is neglected for the sake of simplicity. 
The electronic part is modeled with a first nearest neigh- 
bor single orbital tight-binding (TB) Hamiltonian within 
the orthogonal parametrization [4, 24]. The TB Hamilto- 
nian predicts a zero band gap for all ZGNRs while DFT 
based calculations show that all ZGNRs are semiconduc- 
tors and their band gaps decrease monotonically with 
Nz, for Nz > 4: [4]. For all Nz, highest valence band 
and the lowest conduction band give rise to a high den- 
sity of states near the CNP {E — 0) [■>]. The energy 
band gaps predicted by DFT calculations should not af- 
fect our conclusions because of the transport gap open- 
ing due to disorder. In creating the disordered edges, 
we employ the algorithm explained in Ref. [8]. Edge 
atoms are eroded at each unit cell from both edges. The 
number of carbon atoms to be removed ranges between 
and 4 for each edge. In calculating the transport 
properties, we follow an atomistic approach and employ 
nonequilibrium Green functions within the Landauer for- 
malism [2")], and a decimation technique [::> )] to obtain 
transmission amplitudes Tp\i{u!) and Tc\{E) for phonons 
and electrons, respectively. We refer the reader to Refs. 
[15, 18, 19] for details of the Green function technique. 
We neglect the electron-phonon coupling in this work 
since it is shown that electron-phonon mean free path 
in ZGNRs is tens of /.im at room temperature for ribbons 
having width ^ 10 nm [ ]. The thermoelectric figure 
of merit is defined as ZT = S^GT/k. where S is ther- 
mopower, G is electronic conductance, T is temperature, 
and the thermal conductance with elec- 

tronic and phononic contributions [28, 29]. 
Results and Discussions— We perform electron and 
phonon transport calculations for ZGNRs of two differ- 
ent widths, namely Nz — 10 and 20, with varying the 
length of the disordered region. For each length L, an 
ensemble of edge disordered ZGNRs are generated. Af- 
ter performing the electronic and phononic calculations 
for each configuration, we average the transmission spec- 
tra over the ensembles of 100 disordered edge profiles 




pristine 




■ 


L=15.7nm 






L=63 nm ^ 






L=252 nm 






L=504 nm 

















200 400 600 

Temperature (K) 



200 400 
L (nm) 



600 



FIG. 2: (Color online) The ratio of phonon thermal con- 
ductance fi:ph(r, i) to its pristine value Kp[j(r) is plotted for 
ZGNR(20) at different temperatures. The ratio is fitted to 
a curve using Eq. (1). Inset shows Kph of ZGNR(20) as a 
function of T . 



to obtain Tei{E) and Tph(a'). We quantify the relative 
amount of disorder, Jn^j as the average number of car- 
bon atoms eroded from the edges of a unit cell of GNR 
divided by the number of carbon atoms in a pristine unit 
cell. For the narrow ribbon, the relative amount of dis- 
order is 7io = 0.2, while for the wider one it is 720 = 0.1. 

Ensemble averaged transmission spectra of phonons, 
Tph(ti^), are given in Fig. 1(b) and (c) with varying 
lengths of the disordered region for A'r = 10 and 20, 
respectively. Pristine GNRs display staircase like trans- 
mission spectra, the transmission values corresponding 
to the number of available transport channels, A^ph- The 
transmission values drop significantly even for samples 
as short as 50 nm, the reduction is dramatic for longer 
samples. Comparing the reduction of transmission val- 
ues for two different widths having the same length, one 
observes that higher values result in stronger sup- 
pression of phonon thermal transport. One should note 
that different disorder types can give rise to differences in 
transmission reduction. Our calculations show that edge 
disorder reduces phonon transport effectively for all en- 
ergy values except very low energies. On the other hand, 
isotopic disorder in carbon and BN nanotubes is shown to 
suppress highest energy modes more strongly than other 
modes [12, 13, 14]. Likewise, Anderson-like disorder also 
distinguishes the high energy modes by suppressing them 
more effectively in carbon nanotubes [lo]. Using the re- 
lation Tph(ijj) = Np\i/{1 + L/£p\i) we calculate the elastic 
phonon mean free paths, ^ph. In the insets of Fig. 1, 
£ph are given as functions of energy for Nz = 10 and 20. 
Low energy modes preserve their quasi-ballistic behavior 
through edge disordered zigzag GNRs like they do under 
different disorder types [30]. For a; > 50 cm^^, £ph are 
quite short and always less than 75 nm for both Nz val- 
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ues. Such short mean free paths are observed only at high 
frequencies for isotopic or Anderson-hke disorder. The 
dramatic suppression of phonon transmission can be un- 
derstood in the foUowing way. Anderson-hke or isotopic 
disorder modifies the force constants which gives rise to 
elastic scatterings. Edge disorder, on the other hand, 
not only modifies the vibrational frequencies but also 
changes the number of modes throughout the ribbon. At 
an interface where the width of the ribbon changes, the 
number of transmission channels the ribbon can support 
is changed suddenly (Fig. 1). The difference in trans- 
mission coefficients is significantly large and distributed 
quite homogeneously at all energy values except lu 0. 
This is the reason why €ph is of the same order of mag- 
nitude for all phonon spectrum except at a; 0. Phonon 
thermal conductance, Kph(T). is suppressed strongly with 
increasing L (Fig. 2). Note that Kph is already reduced 
by an order of magnitude for ZGNR(20) at room tem- 
perature when L ~ 250 nm. The length dependence of 
Kph is stronger for the narrow ribbon, i.e. for larger 7. In 
Fig. 2, we plot the ratio of length dependent Kph to its 
pristine value, k°^^, at different temperatures. We observe 
that the ratio follows a single curve for all lengths and 
at temperatures T > 50 K. At low temperatures, high 
frequency phonons are filtered out, and only the low fre- 
quency phonons, which have long £ph, conduct heat. At 
higher temperatures this effect is less pronounced. Since 
£ph is oscillating around similar values except for w 0, a 
definition of an tj— independent effective mean free path, 
£ph, is possible. Letting Tph(t^) ~ iVph(t^)/(l + L/lph), 
one can write 

np,,{T,L)^j^n;,m- (1) 

This approximation accurately reproduces the numerical 
data for L ^ iph, but it is not valid at temperatures T < 
50 K. We perform a fitting and obtain Zph — 23.69 nm for 
ZGNR(20) and 7.04 nm for ZGNR(IO). Note that the fit 
is in very good agreement with the computed data and 
we use the above relation for extrapolating Kph (T, L) for 
L > 500 nm. For L ^ iVph^ph, complete localization 
of phonons in a very wide a range of the spectrum, and 
therefore a glass transition, is expected. 
Charge transport through edge disordered ZGNR shows 
interesting features. In Fig. 3, ensemble averaged trans- 
mission spectra of electrons T^i (E) and calculated mean 
free paths, ici{E), are plotted for = 10 and 20. The 
suppression of Td with increasing L gives rise to differ- 
ent behaviors at different energies. Once edge disorder 
is introduced, the transmission drops significantly at the 
CNP and a transport gap opens [7, (S]. The width of 
the gap is determined by the relative amount of disor- 
der as well as the length of the system, the Opening of 
the transport gap is due to the fact that at the CNP the 
states are fully localized at the edges. As E is increased, 
the charge density is dispersed within the ribbon, and 
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FIG. 3: (Color online) Electron transport through edge dis- 
ordered ZGNRs. Ensemble averaged transmission spectra 
Tel are plotted for pristine ZGNRs and with varying sample 
lengths for Nz = 10 (upper) and A''^ = 20 in (lower). Sample 
lengths are L = 7.9, 31.5, 126, 252 nm (and additionally 1 and 
2 ^m for ZGNR(20)). Elastic mean free paths £ci vs energy 
are extracted from the transmission data. Zero of the energy 
is set to the CNP. 



Tci is very weakly affected by edge disorder. It drops 
abruptly close to the band edges where the density of 
states is singular. Outside the FCP, when N^i > 1, the 
number of elastic scattering channels also increases and 
this results in a large reduction of transmission. One 
should note that id is reduced by approximately an or- 
der of magnitude when the amount of disorder is dou- 
bled. With increasing sample length, a large derivative 
of Tci is observed at the edges of the FCP and close to 
the CNP. Such large derivatives are expected at energies 
where tci{E) changes abruptly, and we exploit this fea- 
ture of edge disordered ZGNR for thermoelectric energy 
conversion. 

To achieve a high ZT , a low thermal conductance to- 
gether with a high Seebeck coefficient (large derivative of 
Tel) and a high charge conductance are required [2(S, 29]. 
In ZGNRs, edge disorder suppresses Kph by few orders 
of magnitude. The derivative of Td is large at the CNP, 
but the transport gap is not sufficiently wide to overcome 
the mutual cancellation of electron and hole contributions 
to 5*. Nevertheless, the cancellation decays at the edges 
of FCP with increasing L as disorder suppresses charge 
transport very strongly out of the FCP. ZT increases 
with increasing length for both ZGNRs, and the peaks 
appear close to the energies where the derivative of Td 
is large. In Fig. 4, ZT is plotted at different tempera- 
tures as a function of the chemical potential and it is 
symmetric with respect to the CNP due to electron-hole 
symmetry. The length of the sample is L = 0.25 nm 
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FIG. 4: (Color online) Thermoelectric figure of merit, ZT , 
versus chemical potential at different temperatures T for 
ZGNR(IO) (left) and ZGNR(20) (right) having lengths L = 
0.25 nm and 4 [ixa, respectively. Zero of [i is the CNP and 
the dotted lines represent the FCP edges. 



(i = 4 [iTo) for iV^ ^ 10 (N^ ^ 20). At room tempera- 
ture ZT reaches the values of 0.39 and 4 for ZGNR(IO) 
and ZGNR(20), respectively. For ZGNR(IO) ZT in- 
creases with temperature up to 1.2 at 750 K, whereas 
for the wider ribbon it decreases with temperature. The 
inverse behaviors with temperature are due to the differ- 
ences in the electronic transmission at energies beyond 
the FCP edges. For ZGNR(IO), Toi is vanishingly smaU 
for S > 1.1 eV, while for ZGNR(20) finite Td remains 
around E = 0.85 eV, which also gives rise to the satellite 
peak of ZT close to this energy. 

In summary, we have shown that phonon thermal con- 
ductance can be suppressed significantly whereas charge 
transport stays intact within the first conduction plateau 
and a large derivative of the electronic transmission func- 
tion is obtained at the edges of the plateau due to edge 
disorder in zigzag graphene nanoribbons. High ZT can 
be achieved depending on the ribbon width and the rel- 
ative amount of disorder. Furthermore, a sharp peak 
in the transmission spectrum, which maximizes ZT [31], 
can be attained due to the fact that the extension of the 
first conduction plateau decreases with increasing ribbon 
width. 
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